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If the dark matter of our galaxy is composed of nuggets of quarks or antiquarks in a colour 
superconducting phase there will be a small but non-zero flux of these objects through the Earth's 
atmosphere. A nugget of quark matter will deposit only a small fraction of its kinetic energy in 
the atmosphere and is unlikely to be detectable. If however the impacting object is composed of 
antiquarks the energy deposited can be quite large and contain a significant charged particle content. 
These relativistic secondary particles will subsequently be deflected by the earth's magnetic field 
resulting in the emission of synchrotron radiation. This work will argue that this radiation should be 
detectable at radio frequencies and that present and proposed experiments are capable of detecting 
such a signal. 



I. INTRODUCTION 

Several recent experiments have made attempts to de- 
tect the radio emission generated by the extensive air 
shower initiated by the impact of an ultrahigh energy 
cosmic ray on the earth's atmosphere [1],[2], [3]. This 
radio wavelength emission is generated by the deflection 
of secondary particles by the earth's magnetic field which 
produces synchrotron radiation. It is the purpose of this 
paper to demonstrate that such experiments are also ca- 
pable of searching for dark matter in the form of heavy 
quark matter nuggets. In the following section I will offer 
a brief overview of the quark nugget dark matter model 
and it's observational consequences for both galactic ob- 
servations and ground based detectors. Section II ex- 
tends the a previous analysis of the air shower induced by 
a quark nugget passing through the earth's atmosphere 
[4] and determines the synchrotron signal generated by 
such an event. Sections III and IV then determine spe- 
cific observable properties of the synchrotron emission. 



A. dark matter from compact composite objects 

The microscopic nature of the galactic dark matter is 
not yet established. While the majority of dark mat- 
ter models introduce a new particle which is fundamen- 
tally weak in its interactions with visible matter, it is 
only the interaction cross section to mass ratio (cr/M) 
that is observationally constrained. As such, a suffi- 
ciently heavy dark matter candidate may avoid obser- 
vational constraints despite a relatively strong interac- 
tion strength. This is possible due to the fact that it is 
only the dark matter mass density which is measured, a 
heavier dark matter candidate will have a lower number 
density and thus smaller flux through any detector. In 
the dark matter model to be discussed here dark matter 
takes the form of heavy compact composite objects with 
nuclear scale densities and composed of the standard light 
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quarks [5] . Several versions of this model have been pro- 
posed dating back to objects such as stranglets [6]. For 
such objects the current observational constraints on the 
dark matter (cr/M) ratio require a mean baryonic charge 
of at least 10 23 . 

Quark nuggets may be formed at the QCD phase tran- 
sition and are stable over the lifetime of the universe. 
Nugget formation allows for the creation of nuggets of 
both matter and antimatter allowing a single process to 
produce both the dark matter and the matter-antimatter 
asymmetry of the universe. This will occur if the produc- 
tion of antinuggets is favoured by a factor of ~ 1.5 over 
the production of matter nuggets leading to dark mat- 
ter which consists of 2 parts matter nuggets to 3 parts 
antimatter nuggets. The excess matter (not bound in 
nuggets) forms the visible matter of the universe in the 
observed five to one matter to dark matter ratio. Orig- 
inally proposed to explain the baryon asymmetry this 
model has been found to have sever observational conse- 
quences for the galactic spectrum. 

• The quark nuggets are surrounded by an "electro- 
sphere" of positrons the outer layers of which are bound 
with energies at typical atomic scales. The annihilation 
of these positrons with the electrons of the interstellar 
medium will result in a positronium decay line (and as- 
sociated three photon continuum) in regions where both 
the visible and dark matter densities are large. In partic- 
ular one should expect a bllkeV line from the galactic 
centre [7] [8]. Such a spectral feature is in fact observed 
and has been studied by the Integral observatory [9]. 

• Positrons closer to the quark matter surface can carry 
energies up to the nuclear scale. If a galactic electron 
is able to penetrate to a sufficiently large depth it will 
no longer produce the characteristic positronium decay 
spectrum but a direct e~e + — > 27 emission spectrum 
[10] . Precisely modeling the transition between these two 
regimes allows for the determination of the strength of 
the MeV scale emissions relative to that of the 511keV 
line [11]. Observations by the Comptel satellite show 
an excess above the background predicted from galac- 
tic sources at the energy and intensity predicted by this 
model [12]. 

• Galactic protons may also annihilate with the an- 
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timatter comprising the quark nuggets. The annihila- 
tion of a proton within a quark nugget will produce 
hadronic jets which cascade down into lighter modes of 
the quark matter. If the energy in one of these jets 
reaches the quark surface it will excite the most weakly 
bound positron states near the surface. These excited 
positrons rapidly loose energy to the strong electric field 
near the quark surface. This process results in the emis- 
sion of Bremsstrahlung photons at x-ray energies [13]. 
Observations by the Chandra observatory indicate an 
excess in x-ray emissions from the galactic centre in pre- 
cisely the energy range predicted [14]. 

• The annihilation of visible matter within the nuggets 
heats them above the background temperature. The 
thermal spectrum from the nuggets may be predicted 
based on the emission properties of the electrosphere 
along with the annihilation rate at various positions 
within the galaxy [15]. The majority of this thermal en- 
ergy is emitted at the eV scale where it is very difficult 
to observe against the galactic background. However the 
emission spectrum will extend down to the microwave 
scale where it may be responsible for the "wmap haze" 

The emission mechanisms involved in describing the 
full nugget spectrum span a very large range of energies 
from the modified thermal spectrum in the microwave up 
to energies associated with nuclear annihilations which 
are observable as gamma rays. However, they make 
no significant contribution to the galactic spectrum at 
the GeV scale or above, an energy range over which the 
FERMI telescope has placed significant constraints on a 
possible dark matter contribution [17], [18], [19],. While 
the uncertainties involved makes an exact determination 
impossible a best fit to the galactic spectrum favours 
quark nuggets with a mean baryon number of roughly 
B ~ 10 25 . 



B. quark matter in the atmosphere 

While the low number density of quark nuggets re- 
quired to explain the observed dark matter mass den- 
sity implies a relatively small flux of nuggets through the 
earth's atmosphere they may still have observational con- 
sequences. Assuming the nuggets to have roughly nuclear 
scale densities observational constraints imply a mini- 
mum mass of a gram and a radius of ~ 10~ 7 m. Given 
the local dark matter mass density of p ~ lGev/cm 3 
and assuming a mean velocity at the galactic scale v g ~ 
200km/s one may estimate the yearly flux of nuggets 
through the atmosphere as, 

this flux is at a level comparable to that of ultra high 
energy cosmic rays near the GZK cuttoff [20] and con- 
sequently may be detectable with large scale cosmic ray 
observatories. However, unlike standard cosmic rays the 



flux of quark nuggets will experience an annual variation 
as the earth moves through the background of galactic 
matter. This annual variation has been used in dark 
matter detection experiments such as Dama [21]. 

In the case of nuggets of quark matter all energy de- 
posited in the atmosphere comes through through colli- 
sional momentum transfer. Relatively little energy is de- 
posited and the nugget continues through the atmosphere 
with virtually no change in velocity. The observational 
prospect for such an event are very low. If however the 
nugget is composed of anti-quarks (as the majority will be 
in the model under consideration) atmospheric molecules 
will annihilate on contact with the quark matter surface 
resulting in the release of substantial amounts of energy 
into the surrounding atmosphere. While the total anni- 
hilation of a quark nugget with B = 10 25 would release 
10 J of energy only a very small fraction of the nugget 
will annihilate in the time it takes to cross the earth's 
surface. As the nuggets carry sufficient momentum to 
traverse the entire atmosphere the limiting factor is the 
amount of atmospheric mass which they encounter. 

This energy will be deposited in the form of thermal 
radiation from the nugget as well as relativistic particles 
and high energy gamma rays emitted in the nuclear inter- 
actions. The prospects for direct detection of secondary 
particles was discussed in [4] . This work focusses instead 
on the prospect of radio frequency detection. Both the 
thermal spectrum and the geo-synchrotron emission from 
the emitted particles will contribute to the radio spec- 
trum, their relative contributions will be approximated 
in section II. 

The total scale of the radio band spectrum is deter- 
mined by the rate at which atmospheric molecules are 
annihilated within the nugget. This rate was estimated 
in [4] where it was shown to increase exponentially with 
the growing atmospheric density until the point at which 
the thermal energy produced by annihilations is sufficient 
to deflect further incoming matter. At this point the 
annihilation rate reaches an equilibrium point and does 
not increase further. The exact value of the equilibrium 
point depends on details of emission from the quark sur- 
face but should occur when the surface temperature is on 
the order of IQkeV . For this temperature scale the anni- 
hilation rate saturates a few kilometers above the earth's 
surface. Above this height the annihilation rate is sim- 
ply determined by atmospheric density and the physical 
cross section of the nugget. 

Tan = a N v N n at (h) h > h eq (2) 
= o N v N n at {h eq ) h < h eq 

Where a at is the physical cross section of the nugget (on 
the order of 10 _10 cm 2 for the nugget mass range consid- 
ered here), and h eq is the height at which the equilibrium 
annihilation rate is reached. 
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II. RADIO FREQUENCY EMISSION 

A. thermal emission 

A quark nugget passing through the atmosphere will 
generate a thermal spectrum which will extend across 
a wide range of energies including the radio band. 
The thermal emission spectrum was first calculated for 
nuggets in the galactic centre [15] where the thermal 
spectrum was calculated to be, 
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This spectrum is modified from a simple black body 
spectrum by the presence of an electrosphere surround- 
ing the quark matter and is suppressed with respect to 
black body radiation for temperatures below the elec- 
tron mass. The nugget temperature is determined by 
the rate at which matter is annihilated within the nugget 
and, as such, will be much higher for a nugget in the 
earth's atmosphere than for nuggets in the interstellar 
medium. The evolution of a nugget's temperature as it 
moves through the atmosphere was described in [4] where 
it was found that the temperature typically peaks near 
lOkeV . At this temperature the emission of radiation in 
the radio band is well described by the w < T limit and 
the spectrum may be taken as 
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In this form it is clear that thermal emission generates 
a relatively flat (log dependance) spectral contribution 
across all radio frequencies. Thermal emission will occur 
uniformly over the nugget surface so that the intensity 
of the thermal component of the spectrum is simply ob- 
tained by dividing expression (4) by Air times the distance 
between the nugget and the observer. 



B. geo-synchrotron emission 

As a quark nugget moves through the atmosphere nu- 
clear annihilations generate a large number of secondary 
particles. Some fraction of these particles, dominated 
by relativistic muons, escape from the nugget into the 
atmosphere. These secondary particles are deflected by 
earth's magnetic field generating synchrotron radiation 
with synchrotron frequency lob = qBo/m. For relativis- 
tic particles in the earth's magnetic field the synchrotron 
radius is much longer than the total path length so that 
the the deflection of the particles is relatively small. In 
this limit we can linearize the equations of motion of a 
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particle within a constant magnetic field. 

v(t) « v + v t « v + ( — ) v x B 

Where vo is the initial velocity of the charged particle 
which, in the case of muons produced in QCD scale pro- 
cesses at the quark surface should be on the order of 
vo ~ 0.9c. The acceleration term in this expression will 
lead to the production of synchrotron radiation. The re- 
sulting radiation is beamed along the direction of the ve- 
locity so that (in the small deflection limit) we need only 
consider the radiation emitted by particles with initial 
velocities directed towards the observer. 

Unlike a conventional air shower the dominant mech- 
anism generating relativistic particles is not direct pair 
production but complex many body interactions at the 
quark matter surface. In the model considered here these 
processes may be summarized by two simple parameters 
the mean velocity of the muons emitted (/?o) and the 
number of muons emitted per nuclear annihilation 
While they could in principle be calculated with in a 
particular quark matter model the introduction of these 
parameters, intended only to capture the approximate 
scale of the radio signal, allows us to discuss the general 
properties of the synchrotron emission without making 
detailed nuclear calculations. Muons produced through 
nuclear interactions within the quark nugget necessarily 
carry nuclear scale energies and thus fio ~ 0.9 — 0.99. The 
muon production coefficient is more difficult to estimate 
but, in the simplified picture where a nuclear annihilation 
in the quark matter produces a muon with properties rel- 
atively close to their vacuum values, it may be estimated 
that only those muons emitted in the direction of the sur- 
face will escape (the remainder being thermalized within 
the nugget) the sharpness of the quark surface also results 
in a high likelihood of an outgoing muon being reflected 
back into the quark matter. The combination of these 
two effects suggests a muon production coefficient on the 
order / M - 0.1 - 0.01. 

While muon production does not occur through direct 
pair production the process is charge independent so fi + 
and fi~ production will proceed at the same rate. In this 
case it is relatively simple to evaluate the electric field 
generated by a charge neutral muon pair in the small lob 
limit. This field at leading order is, 
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Where sin9 V B is the angle between the initial muon ve- 
locity and the earth's magnetic field and R(t) is the dis- 
tance between the muon and the observation point. This 
may be transformed into momentum space as, 
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With Rq being the nugget to observer distance and AR 
the smallest separation between the emitting muon and 
the observer. Note that the relevant timescale in de- 
termining the frequency dependence of the synchrotron 
radiation is the length of time for which the emitting 
particle pair remains relativistic. It should also be noted 
that all information about the shower geometry is carried 
by the sine function and the unitless integral. This ex- 
pression then allows us to estimate the scale of the field 
strength based purely on the numerical coefficient. 
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As should be expected the contribution of a single parti- 
cle pair is relatively weak. In order to determine the total 
field strength associated with the shower this value must 
be scaled up by the total number of particles contributing 
to the radio emission. 

The basic properties of a nugget induced air shower can 
be demonstrated in the geometrically simplified case of a 
shower initiated by quark nugget crossing the atmosphere 
vertically. In this case the earth's magnetic field may be 
taken to lie in the x-z plane with an inclination angle 9b 
and the observation position on the earth's surface by a 
distance from the shower centre (b) and a angle relative 
to the horizontal component of the magnetic field (</>). 
Due to the strong beaming of the synchrotron radiation 
only muons with an initial velocity directed very nearly 
towards the observer need be considered so the particle to 
observer distance may be approximated as R(t) « Ro — 
Poet with i?o = \Jb? + b 2 where h is the nugget's height 
at the time the muon is emitted. This approximation 
breaks down if the muons actual reach the detector, and a 
short distance cutoff must be added for points sufficiently 
close to the shower core. 



III. SIGNAL PROPERTIES 

A. Electric field Magnitude 

The expression (7) gives the momentum space electric 
field strength resulting from a single muon pair. In order 
to determine the total field strength at the earth's surface 
we now sum over all muon pairs emitted towards the 
observer at a given moment. The nuclear annihilation 
rate for a given height is given in expression (2) so that 
the total number of muons produced is determined by 
multiplying this rate by the muon production coefficient 
(/mu) a s discussed above. The rate of muon production 
from a specific point on the nugget surface depends on 
the local flux of material onto the quark surface which 



is proportional to the area perpendicular to the nugget's 
velocity. 
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Here <f> is the angle between the direction of emission and 
the nugget's velocity. As stated above, muon emission 
is dominantly perpendicular to the surface so that the 
angular position of emission on the nugget's surface fully 
determines the direction of the resulting synchrotron ra- 
diation. As such, in the vertical shower case, the radia- 
tion detected a distance b from the shower centre when 
the nugget is a a height hp? has 
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In order to determine the total number of muons con- 
tributing to the observed flux we must estimate the solid 
angle at the nugget's surface which contributes to the 
synchrotron radiation along a given line of sight. To do 
this we note that the intensity of the radiation has an 
angular dependence which scales as S ~ (1 — /3cos0) 4 
which, for /3 close to one, is sharply peaked around zero. 
As such I will take as the angular scale of emission the 
angle at which the intensity falls to half its peak value. 
This angle is defined by the expression 
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Where the second expression uses the small angle approx- 
imation. The relevant solid angle of the nugget surface is 
then given by dfl rs Q\ ; 2 and the total number of muons 
emitted towards the observer at a given time is, 
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Finally we note that the nugget's velocity is much smaller 
than that of the emitted muons so that its height changes 
very little over the time scale on which the muons emit 
synchrotron radiation. This time scale depends on inter- 
action rates with the surrounding atmosphere and a full 
calculation of energy loss rate for a muon is rather com- 
plicated. In order to estimate the time scale involved I 
will simply note that muons lose energy to the surround- 
ing atmosphere more slowly than electrons and assume 
that the muon energy is roughly constant until it decays 
to an electron at which point it is rapidly stopped as 
it scatters off atmospheric molecules. In this approxi- 
mation the muon lifetime sets the timescale over which 
the emitting particle remains relativistic. Under these 
approximation the number of muons contributing to the 
radiation flux when the nugget is at a given height may 
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FIG. 1: Total electric field strength in units of 
(fiV m _1 Mhz~ ) as observed at angle = 250m from 
the shower core when the nugget is at height h=500m (), h = 
lkm (dotted), h = 1500m (dashed). 



FIG. 2: Intensity contours for various magnetic field config- 
urations. The earth's magnetic field is taken to lie in the 
x-z plane with inclination angle <j> = 7r/6, 7t/4, — 7r/6 and 7r/2 
(clockwise from top left.) Distance from shower core is given 
in km. 



be approximated as, 
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where r M is the muon lifetime and 7 is the initial muon 
boost factor. Multiplying this expression by the electric 
field contribution from a single muon pair (7) will then 
give the total field strength when the nugget is at a given 
height. The field strength at a variety of observation 
positions is plotted in figure 1. 



B. geometry 

We now turn to the spacial distribution of the radiation 
as observed at the earth's surface. In the electric field in- 
tensities calculated above the geometry is determined by 
the sin9 V B arising from the cross product of the particle 
velocity with the magnetic field direction. This factor 
favours the emission of synchrotron radiation along lines 
of sight orthogonal to the magnetic field. Along these 
directions the particle deflection, and thus radiation, are 
strongest and the surface radiation acquires a radial de- 
pendence. As the field calculated in (6) resulted from the 
emission of a particle pair with zero net charge the surface 
radiation will be symmetric about the surface projection 
of the earth's magnetic field. The electric field intensity 
at the earth's surface is plotted in 2 for a variety of mag- 
netic field inclinations. 



IV. TOTAL INTENSITY 



The magnitude of the Poynting vector is given by, 
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Or, in momentum space 
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With the field strength for a single muon pair given by 
(7). This radiation is emitted as long as the muon re- 
mains relativistic, as argued above this timescale may be 
roughly estimated by the observer frame muon lifetime. 
In this case the observed flux per charge pair can be ap- 
proximated as, 
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The electric field generated by all particles moving to- 
wards the observer at a given time add coherently so that 
this may be translated to a total flux from all particles by 
multiplying the field strength by the total number of con- 
tributing muon pairs which is half the value determined 
in (14). 
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Finally I want to determine the total intensity resulting 
from the passage of a nugget through the atmosphere. 



I want to add to this the intensity contribution from the 
thermal spectrum (3). As the thermal radiation is emit- 
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frequency (Mhzl 

FIG. 3: Total intensity ( Js~ 1 rn~ 2 Mhz -1 ) received on the 
shower axis (b=0) when the nugget is at heights 1km, 2km 
and 3km (blue, green red respectively.) In each case the ther- 
mal intensity (dashed) is found to be at least two orders of 
magnitude below the contribution from geo-synchrotron radi- 
ation (solid). 
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time (s) 

FIG. 4: Total intensity (Js~ 1 m~ 2 Mhz~ 1 ) received on the 
shower axis (b=0) as a function of time where the nugget is 
taken to reach the surface at t=0. Intensity profiles are shown 
at 5Mhz (blue), 20Mhz(green) and 60Mhz (red). 



ted uniformly from the nugget's surface the total inten- 
sity at a given position is simply given by 



duj dt dA 4?r(6 2 + h 2 N ) du dt y ' 

The intensity of both the thermal and synchrotron ra- 
diation components are plotted in 3. In the range of 
frequencies relevant to most current and planned experi- 
ments the thermal emission has only a minimal impact on 
the total intensity. If however observations are extended 
into the Ghz range the thermal component, with its very 
flat spectrum will become increasingly important. 

Finally, it should be noted that the delivered intensity 
evolves on a timescale much longer than that generated 
by an air shower triggered by a single ultrahigh energy 
proton or nuclei. The timescale in this model is set by 
the time it takes for the quark nugget to pass through the 
earth's atmosphere. Assuming that the quark nuggets 
carry velocities at a typical galactic scale of ~ 200km/ s 
the rise in intensity will take place over tens of millisec- 
onds. Figure 4 shows the rise in intensity as a function 
of time. 



V. CONCLUSION 

This paper has aimed to demonstrate the feasibility 
of detecting the presence of dark matter in the form 
of heavy nuggets of quark matter using radio frequency 
detectors. The passage of a quark nugget through the 
atmosphere will induce an extensive air showers involv- 
ing many secondary charged particles which subsequently 
emit synchrotron radiation across the Mhz band as they 
are deflected by the earth's magnetic field. The resulting 
radio emission is likely to be detectable up to a few kilo- 
meters from the shower core. As such these events should 



be readily detectable by experiments intended to observe 
radio emission from ultra high energy cosmic rays. 

While the intensity of these and event rate of these 
showers may be at a similar scale to that of air showers 
initiated by a single ultra high energy proton or nucleus 
antiquark nugget induced events have several easily ob- 
servable distinguishing properties. 

The nuggets require a time scale on the order of tens of 
milliseconds to traverse the earth's atmosphere and will 
generate a synchrotron signal over much of this time. As 
such the radio signal will be have a duration much longer 
than that of typical cosmic ray events which evolve on 
time scales orders of magnitude faster. 

The nuggets carry galactic scale velocities and their 
flux will show a seasonal variation, however, any micro- 
scopic particle with so small a velocity will be insuffi- 
ciently energetic to initiate an extensive air shower. Con- 
sequently, the detection of a seasonal variation in the air 
shower rate would be a strong indicator of a quark nugget 
contribution to this flux. 

The composite nature of the primary particle in a 
quark nugget initiated air shower means that there is 
a thermal component to the spectrum in addition to the 
geosynchrotron emission generated by the secondary par- 
ticles. While the thermal component is significantly lower 
than the synchrotron signal at the frequencies typically 
observed it's contribution increases at higher frequencies 
and may well be observable as distinct from the syn- 
chrotron signal. 

If, as argued above, large scale cosmic ray detectors 
are also capable of observing the air showers induced by 
dark matter in the form of heavy nuggets of quark matter 
than these properties will allow the two components to 
be readily distinguished through observations at radio 
frequencies. 
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